A multilayered effective medium model for the roughness effect on the Casimir force 
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A multilayered effective medium model is proposed to calculate the contribution of surface rough- 
ness to the Casimir force. In this model the rough layer has its optical properties derived from an 
effective medium approximation, with the rough layer considered as the mixing of voids and solid 
material. The rough layer can be divided into sublayers consisting of different volume fractions of 
voids and solid material as a function of the roughness surface profile. The Casimir force is then cal- 
culated using the generalizations of the Lifshitz theory for multilayered planar systems. Predictions 
of the Casimir force based on the proposed model are compared with those of well known methods 
of calculation, usually restricted to be used with large scale roughness. It is concluded that the 
effect of short scale roughness as predicted by this model is considerably larger than what could be 
expected from the extrapolation of the results obtained by the other methods. 



The calculation of the Casimir force between bodies 
with microscopically rough surfaces has been gaining in- 
creasing attention in the last few years due to the high 
precision experiments that have been carried out [l|. So 
far, only approximate methods of calculation were pro- 
posed, which have different and limited ranges of applica- 
bility. Two are the most commonly used approaches, the 
pairwise summation (PWS) and the proximity force ap- 
proximation (PFA) . Both approaches are assumed to de- 
liver adequate predictions for surfaces having large scale 
roughness, characterized by low spatial frequency (large 
correlation length A), and small surface gradient. How- 
ever, the range of applicability of the PFA was shown to 
extend down to smaller spatial frequencies when com- 
pared to PWS 0. More sophisticated approximated 
methods, based on more general physical concepts, have 
been developed which have a broader and more clearly 
defined range of applicability. For instance, we have the 
methods developed by Maradudin and Mazur [H, and 
more recently by Maia Neto et al. [4]. The method pro- 
posed by Maia Neto et al. is based on the scattering 
approach using a second-order perturbation theory in 
the roughness amplitude, and it can be applied to sur- 
faces with both the average surface separation d and A 
much larger than the rms roughness amplitude a. This 
is qualitatively the same range of applicability of the 
PFA, however, because the perturbative scattering ap- 
proach (PSA) takes into account physical contributions 
not taken into account by the PFA it delivers more re- 
liable results for rough surfaces with smaller spatial fre- 
quencies. More precisely, the following inequalities deter- 
mine the range of applicability of the PWS, a <C d <C A, 
of PFA, a, d < A, and for the PSA it holds that a < A, d. 
These inequalities implie that all these methods can only 
give reliable results when the modulus of the average sur- 
face gradient is much smaller than unity, the restriction 
being less stringent for the PSA 

In most of the experiments performed so far for the 
precise measurement of the Casimir force, the micro- 
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scopic roughness profiles at the surface of the interact- 
ing bodies were such that the use of the PFA or PWS 
was justified Q. However, in a set of recent experi- 
ments [j| the Casimir force between surfaces with large 
amplitude short scale roughness was measured, evidenc- 
ing a large contribution to the force at small separations 
that could not be accounted by the known methods of cal- 
culation. Therefore, an alternative calculational method 
is required. Further motivation for the development of 
a reliable method to calculate the Casimir force between 
surfaces having short scale roughness comes from the po- 
tential relevance of this force in micro- and nanodevices 
operating with small gaps (g, 6- 8] ranging from tens up to 
a few hundreds of nanometers. The short scale roughness 
at the surfaces of such devices results from the limitations 
in the fabrication processes. For instance, for microde- 
vices such as rotary micromotors made from polycrys- 
talline silicon, typically a rms amplitude a — 18 nm Q 
is found. While still in the nanometer range, a smaller 
roughness amplitude can be expected between surfaces 
separated by small gaps due to the more stringent fabri- 
cation requirements. 

In this work we propose a multilayered effective 
medium model (MEMM) that is intended to deliver re- 
liable predictions of the Casimir force between surfaces 
having short scale roughness requiring only a reasonable 
calculational effort. 



I. EFFECTIVE MEDIUM MODEL 

The MEMM is based upon two approximations. The 
first approximation is to consider the rough layer as an 
effective medium whose optical properties result from the 
mixing of voids (vacuum or air) and the material com- 
prising the solid. This effective medium approximation 
(EMA) has been used successfully over decades to model 
the optical properties of rough layers [l(J. The effec- 
tive complex dielectric function e e ff(u)) can, in princi- 
ple, be calculated using one out of various mixing rules 
presented in the literature ffll EH- We choose to use 
the Bruggeman mixing rule, derived from the Brugge- 
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FIG. 1. Schematic drawing showing the rough surfaces and 
an example of a possible effective multilayer model. 



man model |13[ for spherical inclusions, because it cor- 
rectly predicts e e ^ (ui) for any relative fraction of the 
mixed materials, including either sparse or aggregate ran- 
dom configurations that can mimic the structures actu- 
ally found at rough surfaces. In addition, the Bruggeman 
mixing rule has been shown to be in agreement with the 
experimental optical data extracted from light reflected 
from rough surfaces over a wide range of probed wave- 
lengths [Io|, Hi) the agreement being improved by the 
use of multilayer models. 

For the two-phase mixture considered to model the 
rough layer the effective complex dielectric function e e f* 
is obtained from the Bruggeman mixing rule 
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where the input parameters e v and e s correspond to the 
void and solid complex dielectric functions, and / de- 
notes the volume fraction of the solid. For the voids we 
obviously have e v = 1. The optical properties of the solid 
at the rough layer can be approximated by those of the 
bulk, as has been usually done in the calculation based 
on other approaches Q, or some other experimental or 
theoretical dielectric function that is believed to better 
represent the optical properties of the solid at the rough 
layer. 

For a rough surface characterized by the stochastic 
function h{x,y) 1 which denotes the deviation in the z 
direction from the mean value z = 0, the volume frac- 
tion / is generally a function of z. The functional de- 
pendence of / on z is going to determine the effective 
layer thicknesses and average volume fractions for each 
layer used in order to approximate the continuous vari- 
ation of f(z) as depicted schematically in Fig. [TJ Each 
layer has its corresponding e e ^ (w) obtained by solving 
cq. |T]). This is the second approximation we introduce 
into the model, the continuously varying effective dielec- 
tric function e(w, z) is approximated by a z-independent 
function within each layer. This discretization is a well 
known procedure to solve electromagnetic problems in- 
volving continuously variable inhomogeneous media for 
which there is no analytical solution [15j . 



Therefore, in our model we reduce the task of calcu- 
lating the Casimir force between two rough surfaces sep- 
arated by the average gap d to that of calculating the 
force between multilayered systems separated by an ef- 
fective gap / (see fig. Q]). Currently, this calculation can 
only be performed for the case of planar geometry using 
the expressions for the pressure produced by the Casimir 
effect derived by Tomas [l|| and Raabe et al. [l7| , which 
in the more general case of nonzero temperature T can 
be cast in the form [ljj 



P(l) 



E / dqqK J 



-2kI 



-2kI 



(2) 



where the term for 



= must be multiplied by 1/2. 
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In this equation k = \/Cf(i(m)/c 2 + <Z 2 
subara frequencies £ m = 2mi:k^T '/%. r± — r±(q,£) are 
the generalized Fresnel reflection coefficients for s and p 
polarized waves reflecting from the stack of layers above 
(subscript +) or below (subscript — ) the effective gap I. 
In the case of a vacuum gap e(i£) = 1, however, it is 
worth to note that, in general, k is a function of the di- 
electric function calculated over the imaginary frequency 
axes. The Fresnel reflection coefficients can be easily ob- 
tained, for instance, from the set of recurrence relations 
derived in ref . [13] , and result to be functions of the lay- 
ers thicknesses d n ± and their effective dielectric functions 
e n ±(i£)- 

Before we present the results for the Casimir force pre- 
dicted by the MEMM it is worth to discuss its expected 
range of applicability. The relevant parameters for the 
analysis are a — J (h(x, y) 2 ), A, the average surface gra- 
dient (| Vh(x,y) |) at each surface, and the mean gap 
d. For the sake of brevity in the discussion we assume 
all parameters to be at least approximately the same for 
both interacting surfaces, but the discussion could fo- 
cus, for instance, on the surface with the largest a or 
A. We start by noting that the EMA is considered to 
correctly represent the properties of a dielectric medium 
with inclusions in a host material whose dimensions are 
small compared to the wavelength A of the incident elec- 
tromagnetic wave. It is assumed conservatively that the 
largest dimension (height, radius or sides of the surface 
features), l max , must satisfy l max < A/10 [U) E3- Due 
to the symmetric treatment of inclusions and host mate- 
rial the Bruggeman model has the ability to model the 
electrical response of the aggregates (larger structures 
formed above the percolation threshold of the spherical 
inclusions) actually found at a rough surface. When such 
aggregates are present it is the average size of the sur- 
face structures that becomes relevant. Considering the 
symmetric treatment of voids and solid the average lat- 
eral and vertical dimensions of such structures are of the 
order of A and <x, respectively. Now, we have to con- 
sider that the vacuum electromagnetic modes relevant to 
the Casimir effect are those with a wavelength A satisfy- 
ing A ~ d. Therefore, the following inequality must be 



satisfied l max < 0.1A « O.ld, where L 



Max[a, A]. 



3 



This is a conservative limit on l m nx compared to other 
limits found in the literature |18| and the actual limit 
could be less stringent. For this reason the above in- 
equality can still be considered as an adequate criteria 
when the roughness of both surfaces is relevant. However, 
further theoretical investigations should be performed to 
set more precisely the range of validity of the model. 
To conclude, it is worth to note that differing from the 
other methods of calculation, for the proposed model 
there are no upper limits imposed on | Vh(x,y) |. As 
a consequence, the restriction that it C A does not ap- 
ply, demonstrating that the MEMM can be used for short 
scale roughness. 

For the sake of concreteness let us consider the implica- 
tions of the expected range of validity of the model in the 
context of small gap micro and nanodevices, where the 
precise knowledge of the roughness effect on the Casimir 
force can be more relevant. For such devices it is gen- 
erally the case that the Casimir force becomes relevant 
when d < 100 nm In this case, the restriction 

Imax O.lrf implies that A and a at the rough surface 
must be smaller than approximately 10 nm. Considering 
that a is expected to be of the order of a few nanometers 
results that (| Vh(x, y) |) is of the order unity. While well 
suited for the use of the MEMM, this condition can be 
considered out of the range of applicability of the PWS, 
PFA and PSA. 

While no further considerations should be made in the 
case of dielectric materials, for metals the situation is 
more involved. There are two other restrictions. One re- 
sults from the fact that for the validity of the EMA the 
penetration depth 8 of the electromagnetic waves must be 
large compared to the average dimension of the surface 
features. This restriction results from the exponential 
decay of the electromagnetic fields inside lossy inclusions 
such as metals [Tlj . As a consequence, the following in- 
equality must be satisfied l max < o, where S ~ 20 nm for 
good conductors. This is essentially the same restriction 
expected in the context of small gap micro and nanode- 
vices derived above. The other restriction is related to 
the potential effects of spatial dispersion on the deter- 
mination of the Casimir force 0. The effect of spatial 
dispersion in metals can be neglected if the condition 
Vp/u) = VF^/{2nc) -C lmax is satisfied [l9j], where vf 
is the electron Fermi velocity. This condition assures 
that the electrons are going to oscillate restricted to a 
small region of the effective layer were the dielectric func- 
tion does not vary appreciably. As a reference, for gold 
Vf — 1.4 x 10 6 ms _1 , and this restriction implies the ap- 
proximate inequality lQ~ 3 d <C l ma x- This restriction is 
less stringent than the previous ones. 



II. RESULTS 

In what follows we are going to present illustrative re- 
sults for the Casimir force predicted by the MEMM. We 
consider the planar configuration where two rough sur- 
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FIG. 2. Histogram of the function f(z) for a Gaussian height 
distribution. 



TABLE I. Rough layer models. Presented are the number of 
sublayers, layer thickness and volume fraction. 

Model (t of Sublayers Layer thickness / 

1 1 1.2a 0.5 

2 1 1.6a 0.5 

3 2 0.8a/0.8a 0.35/0.65 

4 3 0.4a/1.2a/0.4a 0.2/0.5/0.8 



faces described by the functions h—(x,y) and h + (x,y), 
have an average separation gap d (See fig. [IJ. In the 
MEMM the only required information regarding the 
rough surface are a and the amplitude probability den- 
sity function p(h). We consider the two surfaces having 
a Gaussian height distribution 



p(h) 



2n<r 



(3) 



with the same rms amplitude a. From p(h) the volume 
fraction (material ratio) function f(z) can be calculated 
from the cumulative distribution function as 120] 



p(z')dz' 
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This function varies monotonically from / = 1 (bulk) 
down to / = (vacuum) as z goes from large negative 
up to large positive values. Most of the variation occurs 
within the interval \z\ < 2a as can be inferred from fig. [2 
where the average values of f(z) within intervals of Az = 
0.2u are plotted. 

For the Gaussian surfaces we are considering, with the 
increase in the number of layers, a good convergence was 
already obtained when the rough layer was modeled by 
three layers, resulting in a system comprised of nine lay- 
ers. We have tested several models, comprised of one or 
more layers laying between the vacuum gap and a semis- 
pace having the optical properties of the bulk solid, and 
choose to present the predictions based on four models. 
The models are presented in table HI and are applied for 
both surfaces h~(x,y) and h + (x,y), resulting in a sym- 
metric system with d n+ = d n -, as well as, /„+ = /„-. 
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FIG. 3. £2 as a function of photon energy E for the 
bulk (black/full), and effective medium with / = 0.8 
(black/dashed), / = 0.5 (gray /dashed), and / = 0.2 
(gray/full). The left and right panels are for Si and Au, re- 
spectively. 



In fig. [3] we present illustrative results for the imagi- 
nary part of e e ff(ui) predicted by the Bruggeman model 
and used to calculate e e ^(z£) for each layer. The results 
are for silicon (Si) and gold (Au), and were obtained us- 
ing the experimental data on both ei = Re[e(w)] and 
€2 = Im[e(cj)] as given by Adachi [2l| and Palik [22j], 
respectively. It can be seen that 62 (w) for silicon varies 
smoothly as / increases from values close to zero, charac- 
terizing the prevail of the contribution from the vacuum, 
up to those of the bulk. In the case of gold a more com- 
plex result is evidenced. Below the percolation threshold 
of / = 1/3, ea(<*>) is essentially that expected for an insu- 
lator (see, for instance, the curve for / = 0.2). Above the 
percolation threshold a metallic behavior is observed. It 
is characterized by the divergence of 62(0;) as ui tends to 
zero. It is also worth to note that for / > 0.8 and / < 0.2 

we have the approximate results (w) ~ t\ mlk { UJ ) an d 
1 f (w) ~ 0, for both metals and semiconductors. This 
fact can be used to determine the most adequate thick- 
ness to be considered for the effective medium region 
since the effect on the Casimir force resulting from the 
regions where / > 0.8 and / < 0.2 can be approximated 
by those of the bulk and vacuum, respectively. That is 
the reason why in the models 2, 3 and 4 we have consid- 
ered a roughness layer with thickness ranging from 1.6a 
up to 2cr, encompassing the region that most affects the 
Casimir force. Model 1 was considered for the sake of 
comparison. 

In fig. 31 we present the roughness correction factor 
r\ r (d) to the Casimir force between two semispaces made 
from Si and Au. This correction factor singles out the 
roughness effect from the total force that includes the 
temperature and finite conductivity corrections. The re- 
sults are for a temperature T = 300 K and a value of 
o" = 5 nm was chosen in order to illustrate the effect 
of large amplitude roughness on the Casimir force. We 
investigated the predicted Casimir pressure for Au con- 
sidering e(cj) described by both the plasma and Drudc 
models. The results presented in fig. 2] are those obtained 
using the plasma model, which is expected to provide the 
correct results at finite temperatures 0. rj r (d) calculated 
using the Drude model at T = K differs only slightly 
(by less than approximately 2%) from these results. Due 



to the potential effects of spatial dispersion at large sepa- 
rations, discussed previously, the calculation of rj r (d) for 
Au was restricted to shorter separations. 

In the calculations based on the Drude model we ex- 
tended the experimental data [22| on e(w) to lower fre- 
quencies adopting oj p = 9.0 eV and 7 = 0.035 eV for 
the plasma and relaxation frequencies. Using eq. ([T]) 
the effective e^ 1 F ( w ) IS calculated directly from e^{ui) 
and e 2 9 (w) predicted by the Drude model for the bulk 
and introduced into the usual Kramers-Kronig relation 
giving e e ff(i£,) for the effective layers. For the plasma 
model, the calculations are more involved. In this case 
e|(w) predicted by the plasma model for the bulk is ob- 
tained by subtracting from ef(w) the contribution from 
the relaxation of the conduction electrons. The avail- 
able experimental data for crystalline gold presents in- 
consistencies because they result from the combination 
of different experimental data sets [22[ . For this reason, 
in order to generate a physically consistent e p (w), e\{ui) 
was obtained from e^fw) using the generalized Kramers- 
Kronig relation [23] 



-P 



(5) 



with ujp given above. Finally, e% (w) is determined from 
cq. ((T|) and used to calculate by means of the 

generalized Kramers-Kronig relation whenever / > 1/3. 
This procedure is necessary because of the resulting 
metallic behavior of the effective medium with the associ- 
ated plasma frequency given by {uj^f ) 2 = (3/ — 1)oj^/2. 
For a volume fraction below the percolation threshold the 
effective medium has an insulatorlike dielectric function, 
and the usual Kramers-Kronig relation can be used in 
order to calculate e e ^ 



III. COMPARISON WITH PFA AND PSA 

For the sake of comparison r) r (d) predicted based on 
the PFA and PSA are also presented in fig. 01 While 
approximate analytical expressions for r) PFA (d) can be 
derived, we resorted to the numerical calculation of 



n 



PFA 



(d). A set of synthetic Gaussian surfaces was gen- 
erated and r) PFA (d) calculated for an ensemble of surface 
pairs according to the following expression 

where Sh = 5h{x,y) = h+(x,y) — h-(x,y), A is the sur- 
face area, and rjp{d) is the finite conductivity and tem- 
perature correction factor for the pressure between two 
semispaces having plane surfaces. The use of synthetic 
surfaces allowed us to establish a more realistic scenario 
setting, for instance, the distance of contact between sur- 
faces, which for the ensemble of surface pairs with a = 5 
nm was d ~ 33 nm. For this reason, our results are 



5 



1.8 




50 100 150 200 250 300 




d{nm) 

FIG. 4. r/ r (d) for the models in table [I] model 1 
(black/dot-dashed), 2 (black/dashed), 3 (black/dotted), and 
4 (black/full). Also presented are the predictions based on 
the PFA (gray /full) and PSA (gray /dashed). The upper and 
lower panels are for Si and Au, respectively. 

restricted to d > 35 nm. We note that for stochastic 
rough surfaces the Casimir force and, consequently, r\ r 
predicted by the PFA are only functions of the amplitude 
probability density of the interacting rough surfaces (see 
Section 17.2.2 of ref. @) as for the MEMM. However, 
while these two methods have no explicit dependence on 
other roughness parameters, both A and (| Wh(x,y) |), 
for instance, must be known and taken into account to 
determine which method could provide the most reliable 
predictions. 

From the results presented in fig. 0] we conclude that 
the effects of short scale roughness are underestimated by 
the PFA when compared to the MEMM. Such discrep- 
ancy should be expected due to the distinct ranges of va- 
lidity of these approaches. While in both approaches the 
surfaces are described solely by the same amplitude prob- 
ability density, the actual roughness profiles that could 
be accounted by each model are quite different. For the 
MEMM the surface must have a much smaller A (more 
specifically A <C d) corresponding to a much more com- 
pact roughness profile, which is seem by the relevant elec- 
tromagnetic waves as a flat layer composed of an effective 
material with physical properties between those of the 
vacuum and the bulk. For the validity of PFA, the rough 
surface must be considered as piecewise plane, each piece 
being described as two interacting semispaces, a physical 
picture that differs quite significantly from that consid- 
ered for the MEMM. 



For a comparison with the PSA we resorted to approx- 
imate analytical results. The roughness correction factor 
r]r SA was calculated based on the approximate analytical 
result for the roughness correction predicted for a metal 
described by the plasma model, namely 

^ SA (d)=l + ^A. (7) 

In this equation the Casimir energy relative correction 
factor SE/E = A = 2.7 \f^o 2 j (Ad) was derived in ref. 0] 
under the condition A <C d <C A p , where X p = 136 nm 
is the plasma wavelength of gold. The same condition is 
valid for eq. (J7J, which in the case of a surface covered 
with short scale roughness can only be approximately 
satisfied if the condition a <C A is to be kept. Only as an 
illustrative result, we push the predictions based on the 
PSA slightly beyond the limits of its range of applicability 
and keeping a — 5 nm we assume a short scale roughness 
with A = 3cr = 15 nm in order to plot the curve presented 
in fig. @] Considering the range of validity of eq. (J7J the 
result is approximately valid in the region around d = 60 
nm. For the sake of comparison, while the PSA predicts 
rf(d) oc fitting rf(d) for model 4 leads to the con- 
clusion that rf(d) cx d" 1 - 6 for Au. The exponent a = 1.6 
is between those predicted by the PSA (1 < a < 2) for 
different ranges of the relevant parameters |4{ . Further- 
more, reproducing the trends observed in ref. [1], eq.® 
predicts a roughness correction larger than that of PFA 
in a wide range of separations. In fact, the approximate 
result eq. ((7]) clearly demonstrates that the smaller the 
A the larger the roughness correction, a trend also evi- 
denced in ref. [J by numerical calculations and further 
approximate analytical results. However, for the chosen 
values of a and A the condition (| h(x,y) \) -C 1 can 
not be adequately satisfied and we should rely only on 
the prediction based on the MEMM. Following the trend 
indicated by the approximate PSA results the MEMM 
predicts an even larger roughness correction. This result 
evidences the actual relevance of short scale roughness to 
the Casimir force. 

Due to the importance of correctly establishing the 
range of applicability of the MEMM we further compare 
its predictions with those of the PSA. In order to ob- 
tain reliable predictions from the approximate result eq. 
(|7|) over a wider range of the separation d we consider 
a surface with a smaller roughness amplitude. In fig. [5] 
we present the predictions from the PFA, PSA and the 
MEMM for gold surfaces having a = 2 nm. Therefore, 
the condition A « d « A p can be more appropriately 
satisfied over a wider range of separations d, simultane- 
ously with the constraint a <C A, also required for the 
validity of the PSA. In fig. [5] we present r]* (d) calcu- 
lated for A = 3fi, 4cr and 5(7. For such values of a and 
A the approximate results for the PSA are valid in the 
range 20 < d < 60 nm. The approximate results of the 
PSA indicate a qualitative convergence towards the pre- 
diction of the MEMM before the condition a <C A ceases 
completely to be valid. Therefore, the comparison be- 
tween the predictions of PSA and that of the MEMM 
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FIG. 5. T) r {d) for a Gaussian surface roughness with a — 2 
nm. The curves are the predictions based on the model 4 for 
the MEMM (black/full), PSA with A = 3a (dashed), A = 4a 
(dot-dashed), and A = 5a (dotted), and PFA (gray /full). 



indicates that the proposed model can be used whenever 
A < 3a, well within the range of definition of short scale 
roughness. We can expect that the comparison between 
the predictions based on improved models, such as higher 
order PSA, and the MEMM will give further information 
regarding the range of validity of the proposed model. 



IV. CONCLUSION 

To conclude, it is worth to note that the predictions 
based on the single layer model (model 2) are approxi- 
mately the same as those predicted by the more complex 
three layer model (model 4). This conclusion, being the 
same for Si and Au, suggests that even a single layer 
model can accurately represent the effective properties 
of rough layers. Therefore, with a relatively small calcu- 
lational effort, the MEMM can deliver accurate predic- 
tions of the Casimir force when short scale roughness is 
involved. While the large roughness correction predicted 
by the MEMM is in qualitative agreement with the ex- 
perimental results of ref. [f|, further theoretical and ex- 
perimental investigations are required in order to stablish 
the range of validity of the model. 
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